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A review of coherent phenomena in photoexcited semiconductors is presented. In 
particular, two classes of phenomena are considered: On the one hand the role played 
by optically-induced phase coherence in the ultrafast spectroscopy of semiconductors; 
On the other hand the Coulomb-induced effects on the coherent optical response of 
low-dimensional structures. 

All the phenomena discussed in the paper are analyzed in terms of a theoretical 
framework based on the density-matrix formalism. Due to its generality, this quantum- 
kinetic approach allows a realistic description of coherent as well as incoherent, i.e. 
phase-breaking, processes, thus providing quantitative information on the coupled — 
coherent vs. incoherent — carrier dynamics in photoexcited semiconductors. 

The primary goal of the paper is to discuss the concept of quantum-mechanical phase 
coherence as well as its relevance and implications on semiconductor physics and tech- 
nology. In particular, we will discuss the dominant role played by optically induced 
phase coherence on the process of carrier photogeneration and relaxation in bulk sys- 
tems. We will then review typical field-induced coherent phenomena in semiconductor 
superlattices such as Bloch oscillations and Wannier-Stark localization. Finally, we will 
discuss the dominant role played by Coulomb correlation on the linear and non-linear 
optical spectra of realistic quantum-wire structures. 
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I. INTRODUCTION 



The resonant excitation by coherent optical radiation of an electronic transition in a semiconductor, 
e.g. a valence- to conduction-band excitation, creates a quantum-mechanical coherent superposition of 
the initial and final states of the transition, called optical polarization. The non-linear optical properties 
of this coherent superposition, together with its time evolution, can be used|-|to provide a sensitive 
measurement of many fundamental parameters of the semiconductor materialcl including elastic and 
inelastic scattering, energy level splittings between nearly degenerate states, energy relaxation, as well 
as associated information such as Lande g-factors and band-mixing. 

In general, nonJiaear laser spectroscopy is a very established and well understood means to probe 
these phenomenao However, the description of the non-linear optical response of a semiconduc- 
tor crystal can be considerably more complex than for simple isolated and non- interacting atomsB. 
This is particularly true for the case of the ultrafast optical spectroscopy used for the study of 
the sub-picosecond, jdynamics of photoexcited carriers in bulk systems as well as in semiconductor 
heterostructureadll3'Ell. 

The life-time of the coherent quantum-mechanical superposition generated by an ultrafast laser exci- 
tation, called "dephasing time" , determines the typical time-scale on which coherent phenomena can be 
observed. Such dephasing time reflects the role played by the various "incoherent" , i.e. phase-breaking, 
mechanisms in destroying the phase coherence induced by the laser photoexcitation. Since semicon- 
ductors are chaxacterized by very short electron-hole dephasing times, of the order of few hundreds 
of femtosecondsO, coherent phenomena, manifest themselves only through ultrafast optical experiments 
with sub-picosecond time resolutions!. On this time-scale, the ultrafast evolution of photoexcited 
electron-hole pairs will reflect the strong coupling between coherent and incoherent dynamics, thus 
providing invaluable information on the non-equilibrium relaxation processes occurring in the semicon- 
ductor, e.g. carrier-carrier and carrier-phonon scattering. 

The aim of this paper is to review the basic aspects related to coherent phenomena in semiconductors. 
In particular, we will focus on the ultrafast — coherent vs. incoherent — carrier dynamics as well as on 
Coulomb-correlation effects in photoexcited semiconductors. 

The paper is organized as follows. In the remainder of this section, after a brief historical account of 
coherent experiments in solids, we will try to gain more insight into the concept of coherence by intro- 
ducing a simplified description of the light-matter interaction in terms of a two-level model. In section 
H wc will discuss the theoretical approach commonly used for a realistic description of both coher- 
ent and incoherent phenomena in various semiconductor structures, e.g. bulk systems, semiconductor 
superlattices, quantum wells and wires. Section [II deals with ultrafast carrier photoexcitation and 
relaxation in bulk semiconductors; In particular, we will discuss the dominant role played by coherence 
on the carrier photogeneration process. In section [V we will review and discuss typical field-induced 
phenomena in superlattices, namely Bloch oscillations and Wannier-Stark localization, as well as their 
dephasing dynamics. Sectionals devoted to the analysis of the coherent optical response of quasi-one- 
dimensional systems; More specifically, we will discuss the strong modifications induced by Coulomb 
correlation on the linear and non-linear optical spectra of realistic quantum- wire structures. Finally, in 
section VI we will summarize and draw some conclusions. 



A. Historical background 

Coherent phenomena in-atomic and molecular systems have been investigated for a long timcEI. The 
first spin echo experimentlla was performed in 1950 on protons in a water solution of Fe"'""'""'" ions. Pulses 
in the radio frequency range were generated by means of a gated oscillator with pulse widths between 
20 ^s and a few milliseconds. With these pulses dephasing times of the order of 10 ms have been 
measured. 



In the 1960s echo experiments were brought into the visible rangell3o. A Q-switched ruby laser 
produced pulses of approximately 10 ns duration which were used to observe photon echoes from ruby. 
In this case the dephasing times were of the order of 100 ns. 

As already pointed out, for the observation of any coherent dynamics the pulse width has to be 
shorter than the typical dephasing time. Since in semiconductors electron-hole dephasing times are 
much shorter — they are in the range of a few picoseconds down to some femtoseconds — , coherent 
experiments in semiconductors had to wait until the development of suitable lasers able to generate 
sub-picosecond pulses. 



B. Ultrafast spectroscopy in semiconductors 

The physical phenomena governing the ultrafast carrier dynamics in photoexcited semiconductors can 
be divided into two classes: coherent phenomena, i.e. phenomena related to the quantum-mechanical 
phase coherence induced by the laser photoexcitation, and incoherent phenomena, i.e. phenomena 
induced by the various phase-breaking scattering mechanisms. The above classification in terms of 
coherent and incoherent phenomena is not purely academic; It corresponds to rather different experi- 
mental techniques for the investigation of these two different regimes. 

From an historical point of view, the optical spectroscopy in semiconductors started with the analy- 
sis of relatively slow incoherent phenomena (compared to the electron- hole dephasing time-scale). The 
investigation of nonequilibrium carriers started with the analysis of the incoherent energy-relaxation 
processes in the late 1960s using, cw-lasersEJ. In the 1970s pulse sources for the study of photoex- 
cited carriers became availablalj and this initiated time-resolved studies of the energy-relaxation 
process. Many- experiments based on different techniques have been performedo: band-to-band 
luminescenceED'tj which |»onitors the product of electron and hole distribution functions, band-to- 
acceptor luminescencet3~Ell which, umvides information on the electron distribution functions only, and 
pump-and-probe measurementstJ Ej where the measured differential transmission is proportional to 
the sum of electron and hole distribution functions. 

The theoretical analysis of these relaxation phenomena is commonly based on the semiclassical Boltz- 
mann theory. The Boltzmann equations for both electron pad-hole distribution functions are commonly 
solved by means of semiclassical Monte Carlo simulation£3"LJ. 

In addition to the analysis of incoherent energy-relaxation processes, the ultrafast optical spec- 
troscopy has allowed the investigation of coherent phenomena. As already pointed out, a coherent 
laser field creates, in addition to a non-equilibrium carrier distribution, a coherent polarization. The 
investigation of the coherent dynamics in sepiicanductors started in the 1980s. Different aspects_ham 
been investigatecfl: the optical Stark effectEETE^, the d^ihasing of free carrier£J and excitonstTED, 
quantum beats, related to various tvaes of level splittingsES ij,pCharge oscillations in double-quantum- 
well systpinscjO and superlatticeso, many-particle effectsCjUa, and the emission of coherent THz 
radiationEjEl. 

These experiments cannot be analyzed within the framework of the Boltzmann equation. The reason 
is that this coherent polarization reflects a well-defined phase relation between electrons and holes, which 
is neglected within the semiclassical Boltzmann theory. Any proper description requires a quantum- 
mechanical treatment where, in addition to the distribution functions of electrons and holes, also the 
interband polarization is taken into account as an indepeadent variable. Several approacheS|-iisu« 
been used in the literatureEl: Bogoliubov transformationslEl|— ponequilibrium Green's functionso E3, 
band-edge equatipiisJpased on the real-space density matrixLJ, and the density matrix formalism in 
momentum spaca£3~0. 

During the_.last decade the time resolution has been further improved down to few tens of 
femtosecondsE3. On such extremely short time-scale coherent effects can no longer be neglected and 
are found to play a dominant role also for the case of typical incoherent measurements such as time- 
resolved and time-integrated luminescence. In such conditions, the carrier dynamics is the result of a 
strong interplay between coherent and incoherent phenomena. Therefore, the traditional separation |^ 



^ From an historical point of view, the theoretical approaches commonly used for the interpretation of coherent 
phenomena provide a rather qualitative description of phase-breaking processes in terms of phenomenological 
energy-relaxation and dephasing rates. On the other hand, incoherent phenomena were traditionally inves- 
tigated in terms of semiclassical Monte Carlo simulations, which provide a microscopic description of the 



q 



between coherent and incoherent approaches for the theoretical investigation of photoexcited semicon- 
ductors is no longer valid. What is needed is a comprehensive theoretical framework able to describe 
on the same kinetic level both classes of phenomena as well as their mutual coupling. To this purpose, 
a generalized Monte Cado-jnethod for the analysis of both coherent and incoherent phenomena has 
been recently proposeclE3~Lj.|--Tiie spirit of the method is to combine the advantages of the conven- 
tional Monte Carlo approached in treating the incoherent, i.e. phase-breaking, dynamics with the 
strength of a quantum-kinetic approach in describing coherent phenomenaEl. In particular, the coherent 
contributions are evaluated by means of a direct numerical integration while the incoherent ones are 
"sampled" by means of a conventional Monte Carlo simulation in the three-dimensional k-space. 

This theoretical approach has ben applied successfully to the analysis of vadpus ultrafast optical 
experiments, e.g. four-wave-mixing studies of mEpyJapdy effects in bulk GaAsOO and luminescence 
studies of the hot-carrier photogeneration processO^O. 



C. The meaning of coherence 



In order to clarify the concept of phase coherence, let us consider the simplest|^jDJiysical model for the 
description of light-matter interaction, i.e. an optically driven two-level systemElfl. Within a two-level 
picture, a ground state a with energy Ca and an excited state b with energy et are mutually coupled 
by a driving force, e.g. an external field, and/or by their mutual interaction, e.g. Coulomb correlation. 
The two-level-system Hamiltonian 



is the sum of a free-level contribution 



H = Ho + H' (1) 

Ho = eaalfla + Ebalflb (2) 



and of a coupling term 

H' = C/baa[aa + Uabala^, . (3) 

Here, the usual aecond-quantization picture in terms of creation (a^) and destruction (a) operators has 
been introducecE. The two terms forming the coupling Hamiltonian H' will induce transitions from 
state a to b and vice versa according to the coupling constant Ui^^ = U*^. 

As a starting point, let us consider a single-electron system. In the absence of inter level coupling 
{Uah = 0), we have two stationary states, 

\a{t)) = exp (^-^) at|0) , |b(i)) = exp (^-^) aj, |0) (4) 

(|0) being the vacuum state), corresponding to a single electron in level a or b, respectively. On 
the contrary, in the presence of interlevel coupling, the state of the system is, in general, a linear 
superposition of the non-interacting states in equation : 

m))^Ut)\a{r))+c^{t)\h{t)) , (5) 

whose coefficients obey the following equations of motion: 

d Uab /i(ea-eb)i 



dl'^ ^ ^^"^ n ' 

-Cb = —T-exp Ca . (6) 



dt ih \ Ti 

Again, we see that in the absence of interlevel coupling iU — 0) there is no time variation of the 
coefficients, i.e. if the system is prepared in state |a) or |b) it will remain in such eigenstate. On the 
contrary, the interlevel coupling induces a time variation of the coefficients. 



non-equilibrium scattering dynamics. 



A 



The above two-levekmodel provides the, simplest description of hght-matter interaction in atomic 
and molecular systemsB as well as in solidsfl. The coupling term in equation (|^) has the same structure 
of a field-dipole interaction Hamiltonian and the driving force u can be regarded as a coherent light 
field. We will now consider two limiting cases: the ultrashort- and continuous-excitation regimes. 

An ultrashort optical excitation can be described in terms of a delta-like light pulse: Uab{t) — Ti5{t). 
In this case, the equations of motion can be solved analytically: Due to this excitation, the two-level 
system will undergo an instantaneous transition from its ground state {ca = 1, Cb = 0} to the excited 
state {ca = cosa,Cfa = — isina}, with a = j^. Therefore, after the pulse the system will remain in the 
excited state 



\ii^{t)) = (cosa)|a(i)) + {~isu\a)\h{t)) , (7) 

which is a coherent quantum-mechanical superposition of the two non-interacting states. In addition 
to a finite occupation probability |cbp = siv?a of the excited state, there exists a well defined phase 
coherence between the ground- and the excited-state contributions, i.e. a part from their amplitudes, 
the coefficients Ca and Cb differ in phase by ^. This is what is generally meant by optically induced 
phase coherence. 

Let us now consider a continuous optical excitation resonant with our two-level system: C/ab(i) = 
jj^^iujLt Tff[f^}^ jii^j^ = — Ca- Also for this case, equations (^) can be solved analytically. In particular, 
for the initial condition at t = given by the ground state {ca = 1, Cb = 0}, the solution is given by 

\^/j{t)) = cos{^u;Rt)\a{t)) - ismi^ujRt)\b{t)) , (8) 

where cjr = is the so-called Rabi frequency. Compared to the previous case, the continuous 
excitation gives rise to a periodic population and depopulation of the excited state according to |cbp = 
sin^ {^Lu fit) . This purely coherent phenomenon is known as "Rabi-oscillation regime". As for the 
previous case, the excited state in equation (^) reflects a well-defined optically-induced phase coherence. 

As mentioned above, the interlevel-coupling Hamiltonian in equation (|^) may also describe, in addi- 
tion to an optical excitation, some interlevel coupling due, e.g. to Coulomb interaction. Also in this 
case, in the presence of this interlevel coupling the quantum state of the system will result in a coherent 
superposition of the two stationary states in (^. Again, this coupling will induce a well-defined phase 
coherence between the two levels. 

From the above considerations we see that coherent phenomena can be divided into two basic classes: 
optically-induced and Coulomb-induced phenomena. In order to describe the optical response of an 
electron gas within a semiconductor crystal, the above analysis, based on a single two-level system, 
has to be replaced by a statistical-ensemble approach, i.e. a description based on a collection of 
independent two-level systems. Ajpy statistical ensemble of quantum systems is properly described by 



its single-particle density matrix^ 



(a,VaJ , (9) 



where n denotes a generic set of quantum numbers. Its diagonal elements /„ = p„„ provide the average 
occupation numbers while the non-diagonal terms describe the degree of phase coherence between 
states n and n'. For the case of our two-level system, reduces to a two-by-two matrix: the 

diagonal elements /a = Paa and /b ~ ptt describe, respectively, the average occupation of levels a and 
b, while the non-diagonal terms p — pba and p* = p^b reflect the average degree of phase coherence 
between the ground and the excited state. Starting from the Heisenberg equations of motion for the 
creation and destruction operators, i.e. 

" [4'H] , = K,H] , (10) 

we obtain the following equations of motion for the above density-matrix elements: 

4=i^p+%(A-A)^ (11) 

dt iTi in 

They are known as optical Bloch equationsBJl in analogy with the equations first derived by BlochEl for 
the spin systems. This statistical-ensemble description reduces to the previous one for the case of a so 



called "pure state" , i.e. the case in which all the two-level systems forming the ensemble are in the same 
quantum-mechanical state \^). In this case, the above set of kinetic equations and the equations of 
motion (||) are totally equivalent. However, as we will discuss in section |l[ the density-matrix approach 
introduced so far allows, in addition to the study of coherent phenomena, the analysis of incoherent 
phenomena, which is not possible within a simple Schrodinger-equation formalism. 

The optical Bloch equations (^0|) provide the simplest description of light-matter interaction. Let us 
consider again the case of a continuous optical excitation resonant with our two-level system: Uab{t) = 
C/oB*"^*. If we choose as initial condition at time t = the state {/a — /o, ft — 0,P — 0}, the solution 
of the above optical Bloch equations is given by: 

fa{t) = foCOS^i^UJRt) 

fb{t) = foSin^i^UJRt) 



fo 



p{t) ^ ^ [e— ' - e-^'^ 'J , (12) 

where uj^ = ujl i i^r- As for the case of a single two-level system, the above solution describes 
a Rabi-oscillation regime. In particular, the interlevel density-matrix element p originates from the 
superposition of the two frequency components w+ and lu^ . They differ from ujl by the Rabi frequency 
u>R. Such modification of the two- level frequency uj^ due to its coupling with the external field is known 
as "Rabi splitting"El. 

If we now rewrite the interlevel density-matrix element p in ( p^ ) as 

p(t) = -^e— *sin(c.Kt) , (13) 

we see that, a part from the quantum-mechanical phase factor corresponding to the energy separation 
hujL = Efa — Eai its amplitude exhibits Rabi oscillations according to sin(tt'iji). More specifically, we 
obtain: 

\pf (X sin^(w_Rt) oc cos^i^ujRt) sin^(icjflt) oc /a(i)/b(*) . (14) 

i.e. the quantity |pp is proportional to the product of the two occupation numbers, thus rcficcting the 
total (or macroscopic) dipole moment of the two-level system at time t. This elucidates the link between 
optically-induced phase coherence and polarization: a coherent optical excitation gives rise to a coherent 
quantum-mechanical superposition of the two states which results in a macroscopic polarization of the 
system. Such polarization field is fully described by the non-diagonal matrix element p in (p^). 

The above simplified description of light-matter interaction neglects any kind of incoherent, i.e. 
phase-breaking, phenomena. As we will see in section |l^, such dephasing processes will lead to a decay 
of the interband polarization p, thus destroying the optically induced phase coherence in the carrier 
system. 

As a final remark, let us discuss the concept of phase coherence in connection with the choice of 
representation. On the basis of the density-matrix formalism considered so far, the phase coherence 
is described by the non-diagonal density-matrix elements pnn' in (^). However, this separation in 
diagonal and non- diagonal terms is clearly representation dependent: what is diagonal in a given basis 
is in general not diagonal in a different basis and vice versa. If one considers, as basis set, the eigenstates 
of the total Hamiltonian H (which includes the external driving force), the density matrix p in (||) is 
always diagonal, i.e. no phase coherence. Thus, in order to speak of phase coherence, we need to regard 
the total Hamiltonian H as the sum of the system Hamiltonian Ho plus a driving-force term H' (see 
equation (|l|)). Provided such separation between system of interest and driving force, the non-diagonal 
density-matrix elements within the representation given by the eigenstates of Ho will describe the 
degree of coherence induced in the system by the driving force H'. 

Two typical cases, discussed respectively in sections [V and may help in clarifying the meaning 



of coherence. The first case is that of a semiconductor superlattice in the presence of a constant 
and homogeneous electric field: As we will see, within the Wannier-Stark representation the Bloch- 
oscillation dynamics is the result of phase coherence between the various Wannier-Stark states; On 
the contrary, within the Bloch representation the same phenomenon is purely described in terms of 
carrier populations, i.e. no off-diagonal density-matrix elements. The second case is that of Coulomb- 
induced phase coherence: Within a free-particle representation, even at the simplest Hartree-Fock level. 
Coulomb interaction induces phase coherence between free-carrier states. This is clearly not the case 
within an exciton basis, where the density matrix is diagonal. 



II. THEORETICAL BACKGROUND 



In this section we will review in a systematic way the basic ideas used in the theoretical analysis of 
ultrafast carrier dynamics in semiconductors. The approach we are going to present is based on the 
density-matrix formalism introduced in section |[ Cj . 



A. Physical system 

In order to study the optical and transport properties of semiconductor bulk and heterostructures, 
let us consider a gas of carriers in a crystal under the action of an applied electromagnetic field. The 
carriers will experience their mutual interaction as well as the interaction with the phonon modes of 
the crystal. Such physical system can be described by the following Hamiltonian: 



H — He + Hn + He 



H 



cp 



H 



pp 



(15) 



The first term describes the noninteracting-carrier system in the presence of the external electromagnetic 
field while the second one refers to the free-phonon system. The last three terms describe many- 
body contributions: they refer, respectively, to carrier-carrier, carrier-phonon, and phonon-phonon 
interactions. 

In order to discuss their explicit form, let us introduce the usual second-quantization field operators 
4't(r) and '4'(r). They describe, respectively, the creation and the destruction of a carrier in r. In 
terms of the above field operators the carrier Hamiltonian He can be written as 



He y" rfr*^(r) 



(-iWr-fA(r,t))^ 

2mo 



*(r) 



(16) 



Here, F'(r) denotes the periodic potential due to the perfect crystal while A(r, t) and (p{r, t) denote, re- 
spectively, the vector and scalar potentials corresponding to the external electromagnetic field. Since we 
are interested in the electrooptical properties as well as in the ultrafast dynamics of photoexcited carri- 
ers, the electromagnetic field acting on the crystal — and the corresponding electromagnetic potentials — 
will be the sum of two different contributions: the high-frequency laser field responsible for the ultra- 
fast optical excitation and the additional electromagnetic field acting on the photoexcited carriers on 
a longer time-scale. More specifically, by denoting with the labels 1 and 2 these two contributions, we 
can write 



A(r,i) = Ai(r,i) + A2(r,t) , ip{r,t) = <^i(r, i) -t- <^2(r, i) 



and recalling that 



E(r,i) 



ld_ 

cdt 



A(r, t) - Vr¥'(r, t) , B(r, t)=\7,x A(r, t) 



we have 



E(r,i) =Ei(r,t)+E2(r,t) , B(r, t) = Bi(r, t) + B2(r, t) . 



(17) 



(18) 



(19) 



Equation ( p^ ) , which gives the electromagnetic fields in terms of the corresponding vector and scalar 
potentials, reflects the well known gauge freedom: there is an infinite number of possible combinations 
of A and ip which give rise to the same electromagnetic field {E,B}. We will use such freedom of 
choice for the laser field (term 1): we assume a homogeneous (space-independent) laser field Ei(t) fully 
described by the scalar potential 



(^i(r,t) = -Ei(t)-r . 



(20) 



This assumption, which corresponds to the well known dipole approximation, is well justified as long 
as the space-scale of interest is small compared to the light wavelength. The explicit form of the laser 
field considered here is 



E:(t)e 



(21) 



-7 



where Eo{t) is the amphtude of the hght field and ojl denotes its central frequency. 

With this particular choice of the electromagnetic potentials describing the laser field, the Hamilto- 
nian in (|l^) can be rewritten as 



H, = H° + H,, , (22) 



where 



(-zWr-fA2(r,t))' 



+ e(/:>2(r,t) + F'(r) 



*(r) (23) 



2mo 

describes the carrier system in the crystal under the action of the electromagnetic field 2 only, while 

H,,=ey'dr*t(r)^i(r,i)*(r) (24) 

describes the carrier-light (cl) interaction due to the laser photoexcitation. 

In analogy with the carrier system, by denoting with h]^ ^ and ^ the creation and destruction 
operators for a phonon of mode A and wavevector q, the free-phonon Hamiltonian takes the form 

U^p = J2hu;^^bl^b^^, (25) 

qA 

where ujqx is the dispersion relation for the phonon mode A. 

Let us now discuss the explicit form of the many-body contributions. The carrier-carrier interaction 
is described by the two-body Hamiltonian 



^ J dr J dr'*t(r)*t(r')Vcc(r-r')*(r')*(r) , (26) 



where Vcc denotes the Coulomb potential. 

Let us now introduce the carrier-phonon interaction Hamiltonian 



Hep = J dr*t(r)v;p(r)*(r) , (27) 



where 

E 

qA 



[g^xbqxe''''' + r^xbl.e-^^-] (28) 



is the electrostatic phonon potential induced by the lattice vibrations. Here, the explicit form of the 
coupling function depends on the particular phonon mode A (acoustic, optical, etc.) as well as on 
the couphng mechanism considered (deformation potential, polar couphng, etc.). 

Let us finally discuss the phonon-phonon contribution Hpp. The free-phonon Hamiltonian Hp in 
(p5|), which describes a system of noninteracting phonons, by definition accounts only for the harmonic 
part of the lattice potential. However, non-harmonic contributions of the interatomic potential can 
play an important role in determining the lattice dynamics in highly excited systemsEZI, since they are 
responsible for the decay of optical phonons into phonons of lower frequency. In our second-quantization 
picture, these non-harmonic contributions can be described in terms of a phonon-phonon interaction 
which induces, in general, transitions between free-phonon states. Here, we will not discuss the explicit 
form of the phonon-phonon Hamiltonian Hpp responsible for such a decay. We will simply assume that 
such phonon-phonon interaction is so efficient to paaintain the phonon system in thermal equilibrium. 
This corresponds to neglecting hot-phonon effectala. 

It is well known that the coordinate representation used so far is not the most convenient one in 
describing the electron dynamics within a periodic crystal. In general, it is more convenient to employ 
the representation given by the eigenstates of the noninteracting-carrier Hamiltonian — or a part of 
it — since it automatically accounts for some of the symmetries of the system. For the moment we will 
simply consider an orthonormal basis set {(/)„(r)} without specifying which part of the Hamiltonian 
is diagonal in such representation. This will allow us to write down equations valid in any quantum- 
mechanical representation. Since the noninteracting-carrier Hamiltonian is, in general, a function of 



time, also the basis functions (f>n niay be time-dependent. Here, the label n denotes, in general, a set 
of discrete and/or continuous quantum numbers. In the absence of electromagnetic field, the above 
wavefunctions will correspond to the well known Bloch states of the crystal and the index n will reduce 
to the wavevector k plus the band (or subband) index v. In the presence of a homogeneous magnetic 
field the eigenfunctions 0„ may instead correspond to Landau states. Finally, for the case of a constant 
and homogeneous electric field, there exist two equivalent representations: the accelerated Bloch states 
and the Wannier-Stark picture. Such equivalence results to be of crucial importance in understanding 
the relationship between Bloch oscillations and Wannier-Stark localization and, for this reason, it will 



be discussed in more detail in section [V A 



Let us now reconsider the system Hamiltonian introduced so far in terms of such 0„ representa- 
tion. As a starting point, we may expand the second-quantization field operators in terms of the new 
wavefunctions: 



*(r) = ^0„(r)a„ , ^^(r) = ^C(r)at 



(29) 



The above expansion defines the new set of second-quantization operators aj^ and a„; They describe, 
respectively, the creation and destruction of a carrier in state n. 

For the case of a semiconductor structure (the only one considered here) , the energy spectrum of the 
noninteracting-carrier Hamiltonian ( p3|) — or a part of it — is always characterized by two well-separated 
energy regions: the valence and the conduction band. Also in the presence of an applied electromagnetic 
field, the periodic lattice potential V'' gives rise to a large energy gap. Therefore, we deal with two 
energetically well-separated regions, which suggests the introduction of the usual electron-hole picture. 
This corresponds to a separation of the set of states {(f>n} into conduction states {(p^} and valence 
states {0j }■ Thus, also the creation (destruction) operators ajj (a„) introduced in equation ( p9[ ) will 
be divided into creation (destruction) electron and hole operators: 
new electron- hole picture, the expansion (p9|) is given by: 



(Cj) and (d ). In terms of the 



i 3 

*t(r)=^C(r)cl+^0^'(r)d,- . 



(30) 



If we now insert the above expansion into equation ( p3| ) , the noninteracting-carrier Hamiltonian takes 
the form 



n' 



(31) 



where 



2rno 



■ eip2 + V~ Co 



(32) 



are just the matrix elements of the Hamiltonian in the ^-representation. The ± sign refers, respectively, 
to electrons and holes while €0 denotes the conduction-band edge. Here, we neglect any valence-to- 
conduction band coupling due to the external electromagnetic field and vice versa. This is well fulfilled 
for the systems and field-regimes we are going to discuss in this paper. As already pointed out, the 
above Hamiltonian may be time-dependent. We will discuss this aspect in the following section, where 
we will derive our set of kinetic equations. 

Let us now write in terms of our electron-hole representation the carrier-light interaction Hamiltonian 
Ml): 



^4E 



-{MA 



(33) 



The above expression has been obtained within the well known rotating-wave approximatiorn by ne- 
glecting intraband transitions, absent for the case of optical excitations. Here, /Ltfj* denotes the optical 
dipole matrix element between states (p'i and 



Q 



Similarly, the carrier-carrier Hamiltonian (^6|) can be rewritten as 

H, 



^ , „t 



r, ''ill 

ililizii 

E ^^.mA/n^nCi., (34) 



where 

= / '^^ / (r')V^^^(r - ^')4>i. {r')4>u (r) (35) 

are the Coulomb matrix elements within our ^-representation. The first two terms describe the repulsive 
electron-electron and hole-hole interaction while the last one describes the attractive electron-hole 
interaction. Here, we neglect terms that do not-conserve the number of electron-hole pairs, i.e. impact- 
ionization and Auger-recombination processealj, as well as the interband exchange interaction. This 
monopole-monopole approximation is justified as long as the exciton binding energy is small compared 
to the energy gap. 

Finally, let us rewrite the carrier-phonon interaction Hamiltonian introduced in equation (|27|): 



with 



ii' ,qA 

E [4/,qA4V^/+5i;',qA4'<Arf,] (36) 



ml\,^9^J dv^l'^*{r)e^-4>l!\v) . (37) 



In equation (|3^) we can clearly recognize four different contributions corresponding to electron and hole 
phonon absorption and emission. 



B. Kinetic description 

Our kinetic description of the ultrafast carrier dynamics in semiconductors is based on the density- 



/S' = (ck'), /■?,■' = (rfk') (38) 



matrix formalism, introduced in aeciion |I C| . Since this approach has been reviewed and discussed in 
several papersEflEj and text-booksEllj, here we will simply recall in our notation the kinetic equations rel- 
evant for the_analysis of carrier dynamics in semiconductor heterostructures, generalizing the approach 
presented ino to the case of a time-dependent quantum-raechanical representation. 

The set of kinetic variables is the same considered inEJ. Given our electron-hole representation 
{0^}, {0^}, we will consider the intraband electron and hole single-particle density matrices 
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as well as the corresponding interband density matrix 

Pji = {djC,) . (39) 

Here, the diagonal elements ff^ and /j^ correspond to the electron and hole distribution functions of 
the Boltzmann theory while the non-diagonal terms describe intraband polarizations. On the contrary, 
the interband density-matrix elements p^^ describe interband (or optical) polarizations. 

In order to derive the set of kinetic equations, i.e. the equations of motion for the above kinetic 
variables, the standard procedure starts by deriving the equations of motion for the electron and hole 
operators introduced in (|30|): 



in 



= y dr</.r(r)*(r) , = J dr^f (r)*t(r) . 

By applying the Heisenberg equation of motion for the field operator i.e. 



dt ih^ ' ^ 



it is easy to obtain the following equations of motion: 



d 1 r n 

di'^ = ^n^-^ 



dt 



(40) 



(41) 



(42) 



with 



Zf/'' ^ih I dr 



dt 



(43) 



As for the case of equation (|3l|), here we neglect again valence-to-conduction band coupling and vice 
versa. Compared to the more conventional Heisenberg equations of motion, they contain an extra-term, 
the last one. It accounts for the possible time dependence of our ^-representation which will induce 
transitions between different states according to the matrix elements Z^, . 

By combining the above equations of motion with the definitions of the kinetic variables in ( ^8[]39| ) , 
the resulting set of kinetic equations can be schematically written as: 
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(44) 



They exhibit the same structure of the equations of motion ( ^^ for the electron and hole creation and 
destruction operators: a first term induced by the system Hamiltonian H (which does not account for 
the time variation of the basis states) and a second one induced by the time dependence of the basis 
functions (p. 

Let us start discussing this second term, whose explicit form is: 



dt 
d 



■I 1-L12 

df'^'^^ 

d_ 

di^^''' 



in ^ 



1224 ^"^1^3 



J I: 



iTi ^-^ 



1J3 



i2j2 



(45) 



As we will see in section [VA, these contributions play a central role for the description of Zener 
tunneling within the vector-potential representation. 

Let us now come to the first term. This, in turn, is the sum of different contributions, corresponding 
to the various parts of the Hamiltonian. The total Hamiltonian can be regarded as the sum of two 
terms, a single-particle contribution plus a many-body one: 



H = H, 



H,i + Hp) + (H, 



(46) 
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The explicit form of the time evolution due to the single-particle Hamiltonian Hsp (non-interacting 
carriers plus carrier-light interaction plus free phonons) is given by: 




with U,,j,^-fil^^^E-{t). 

This is a closed set of equations, which is a consequence of the single- particle nature of H^p. In 
addition, we stress that the structure of the two contributions entering equation ( ^ ) is very similar: 
one can include the contribution (^) into (^7|) by replacing e with e + Z. 

Let us finally discuss the contributions due to the many-body part of the Hamiltonian: carrier- 
carrier and [Carrier-phonon interactions (the phonon-phonon one is not explicitly considered here). As 
discussed iria, for both interaction mechanisms one can derive a hierarchy of equations involving higher- 
order density matrices. In order to close such equations with respect to our set of kinetic variables, 
approximations are needed. The lowest-order contributions to our equations of motion are given by 
first-order terms in the many-body Hamiltonian: Hartree-Fock level. Since we will neglect coherent- 
phonon states, the only Hartree-Fock contributions will come from carrier-carrier interaction. They 
simply result in a renormalization 

Kt--E^:Uu//t (48) 

of the single-particle energy matrices e*^^^ as well as in a renormalization 

AC/.,,, (49) 

of the external field U. ^ We stress that the Hartree-Fock approximation, which consists in factorizing 
average values of four-point operators into products of two density matrices, is independent from the 
quantum-mechanical picture. It is then clear that the above kinetic equations are valid in any quantum- 
mechanical representation. 

All the contributions to the system dynamics discussed so far describe a fully coherent dynamics, 
i.e. no scattering processes. In order to treat incoherent phenomena, e.g. energy relaxation and 
dephasing, one has to go one step further in the perturbation expansion taking into account also 
second-order contributions (in tiie perturbation Hamiltonian Umb)- The derivation of these higher- 
order contributions, discussed ina, will not be repeated here. Again, in order to obtain a closed set of 
equations (with respect to our set of kinetic variables (|3^p9|)) additional approximations are needed, 
namely the mean-field and Markov approximations. As for the Hartree-Fock case, the mean-field 



^The explicit form of the renormalization terms considered in this paper accounts for the Fock contributions 
only, i.e. no Hartree terms. The general structut* of Hartree-Fock contributions, relevant for the case of a 
strongly non-homogeneous system, is discussed inEj. 
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approximation allows to write the various higher-order density matrices as products of single-particle 
ones. The Markov approximation allows to eliminate the additional higher-order kinetic variables, e.g. 
phonofua^isted density matrices, providing a closed set of equations still local in time, i.e. no memo» 
efFectsH~E3. This last approximation is not performed in the quantum-kin^ic theory discussed inp 
where, in addition to our single-particle variables, one considers two-particleL^ and phonon-assistccEi 
density matrices. 

While the mean-field approximation is representation-independent, this is unfortunately not the case 
for the Markov limit. This clearly implies that the validity of the Markov approximation is strictly 
related to the quantum-mechanical representation considered. 

The above kinetic description, based on intra- and interband density matrices, allows us to evaluate 
any single-particle quantity. In particular, for the analysis of the ultrafast carrier dynamics in pho- 
toexcited semiconductors two physical quantities play a central role: the intra- and interband total (or 
macroscopic) polarizations: 

P<^/\t) = J2 M^^ni^it) . P^'^ = J2 A^?>,.(i) . (50) 

II' V 

where M'^/'' and ^'^^ denote, respectively, the intra- and interband dipole matrix elements in our ^-repre- 
sentation. The time derivative of the intraband polarization P«^/'' describes the radiation field induced 
by the Bloch-oscillation dynamics (which for a superlattice structure is in the TeraHertz range) while 
the Fourier transform of the interband (or optical) polarization P'^'' provides the linear and non-linear 
optical response of the system. 



III. COHERENT ULTRAFAST DYNAMICS IN BULK SEMICONDUCTORS 



In this section, we will discuss the dominant role played by optically-induced phase coherence on 
the ultrafast generation and relaxation of photoexcited carriers in bulk semiconductors. In particular, 
we will compare the description based on the theoretical approach discussed above with the more 
conventional picture based on the Boltzmann theory. 



A. Bloch model 



In order to study the ultrafast carrier dynamics in bulk systems, let us consider a two-band semi- 
conductor model. For the case of a homogeneous system, the only relevant terms of the single-particle 
density matrix in k-space are the diagonal ones. This property, due to the translational symmetry of the 
Hamiltonian, reduces the set of kinetic variables (|3^p9|) to the following electron and hole distribution 
functions (intraband density-matrix elements) 

/k = (4ck) , /-k = {d\d_^) (51) 

together with the corresponding polarizations (interband density-matrix elements) ^ 

Pk = ('^-kCk) ■ (52) 
The explicit form of the kinetic equations ( |4^ ) within the above two-band picture is given by 

5k(t) - E [^k'k/k (1 - f^) - W'kk'/k' (1 - /k)] 

k' 

.9_k(t) - E [<k/k (1 - /k') - <k'/k' (1 - /k)] 

k' 



d 
It 



-ft- 



Here, the standard electron-hole picture introduced in section [I A has been applied to our plane-wave states. 
In particular, due to the charge-conjugation symmetry, the hole states are still labeled in terms of the corre- 
sponding valence-electron states, i.e. k'' = — k'^. 
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J^P^. = ^ [(4 + + f/k (1 - - f\)] 

-Et^'k^'k-^k'Pk'] , (53) 

k' 

with the generation rate 

5k = ^ [C/kPk - f/kPk] ■ (54) 

Here, the incoherent contributions are treated within the usual Markov Hmit as described i 
Within such approximation scheme, the incoherent contributions to the polarization dynamics exhibit 
the same structure as for the distribution functions in terms of the following in- and out-scattering 
rates 

W^k'k = ^ E [^k'k (1 - /k') + ^k'k/k'] . (55) 

v—e^h 

e/h 

where W^,^ are the usual scattering rates of the semiclassical Boltzmann theory. This Boltzmann-like 



structure of the scattering term in the paLaxization equation is the starting point of the generalized 
Monte Carlo method discussed in section |p3T3. 

The above kinetic equations are known as semiconductor Block equations (SBE)^. As we can see, here 
the carrier photogeneration is a two-step process: The external field U induces a coherent polarization 
p which, in turn, generates electron- hole pairs via its coupling with the field according to equation (Js^). 
The generation rate is thus determined by the interband polarization which, in turn, is influenced 
by the various density-dependent scattering mechanisms, e.g. carrier-carrier processes. Therefore, in 
contrast to the semiclassical case discussed below, the generation rate within the SEE model is clearly 
density dependent. 



B. Boltzmann model 

The conventional Boltzmann model is obtained from the above SEE by performing an adiabatic 
elimination of the interband polarization pk, i-e. by inserting a formal solution of the polarization 
equation into the generation rate, ^3) and then performing a Markov limit with respect to the electron 
and hole distribution functionsEj'tJ. This corresponds to assuming that the carrier distributions are 
slowly-varying functions on the time-scale of the laser photoexcitation, i.e. the carrier generation 
and relaxation are treated as independent processes. Within such approximation scheme, except for 
phase-space-filling effects, the generation rate entering the semiclassical Eoltzmann equations (EE) for 
electrons and holes is fully determined by the temporal and spectral characteristics of the laser pulse. 
Contrary to the SEE case, all effects related to the optically induced phase coherence — and to its 
temporal decay — are neglected. 



C. Coherent carrier photogeneration 

In order to compare the differentr«iodels of carrier dynamics in photoexcited semiconductors, 
review some simulated experimentalfl based on the generalized Monte Carlo solution of the SBI 
discussed in section |. The corresnonding BE has been solved, for comparison, by a standard ensemble 
Monte Carlo (EMC) simulationBS. Both BE and SEE have been solved numerically for the case of 
GaAs bulk excited by a 150 fs laser pulse. 

In figure |l| the generation rate at different times as obtained from the SEE is plotted as a function of 
the carrier wave vector for three different densities. At the lowest density the behavior is essentially the 
same as in the case without carrier-carrier scattering: Energy-time uncertainty leads to an initially very 



In the absence of C oulo mb correlation and scattering, the SBE in (B3) reduce to the optical Bloch equations 



introduced in section IC, i.e. they describe a collection of independent (non-interacting) two- level systems. 



broad generation rate; with increasing time the hne narrows and, in the tails, ealiibits negative parts 
due to a stimulated recombination of carriers initially generated ofT-resonanca£3'Ej. After the pulse, 
the distribution function of the generated carriers is in good agreement with the BE result. Scattering 
processes destroy the coherence between electrons and holes which is necessary for the stimulated 
recombination processes. As a consequence, with increasing density the negative tails are strongly 
reduced and the generation remains broad for all times resulting in a much broader carrier distribution 
than in the BE case. 



D. Comparison with experiments 

Hot carrier luminescence has provep-jto be a powerful technique to study the ultrafast dynam ics , o i j P |lji% 
toexcited carriers in semiconductorstS. In band-to-acceptor (BA) luminescence experimentstSrcilO'Qj, 
due to their high sensitivity, carrier densities as low as several 10^'^ cm~^ have been reached. Thus, 
the transition between a dynamics dominated by carrier-phonon scattering to a dynamics dominated 
by carrier-carrier scattering (with its consequences for the generation process discussed above) is ex- 
perimentally accessible. _ _ 

In the BA luminescence experiment reviewed hereLjT3, a 3 /im thick GaAs layer doped with Be 
acceptors of a concentration of 3 x 10^^ cm~^ is excited at a photon energy of 1.73 eV by transform 
limited L5p fs pulses from a mode-locked Ti:sapphire laser. Time-integrated luminescence spectra are 
recorded^. The right column of figure ^ shows measured BA luminescence spectra for three different 
values of the carrier density. At low density we observe an initial peak of the generated carriers 
(marked bold) and pronounced replicas due to the emission of an integer number of optical phonons. 
With increasing density the peaks become broader and at the highest density only a slight structure 
related to the phonons is still visible. 

These measured spectra are compared with corresponding simulated experiments based on the BE and 
SBE modelsLj (left and middle column in figure H) . We find several pronounced differences between the 
two models. In the semiclassical model the unrelaxed peak is clearly visible up to the highest density 
in contrast to the coherent model, where, in agreement with the experiment, this peak is strongly 
broadened. At the higher densities the semiclassical spectra exhibit an increase in the broadening of 
subsequent replicas which is not present in the Bloch-model calculations as well as in the experimental 
results. 

To illustrate this difference quantitatively, in figure ^ we have plotted the full width at half maximum 
(FWHM) of the three highest peaks in the spectra. From the semiclassical model we obtain a strong 
increase in the density dependence of subsequent peaks. This behavior can be easily understood: Due 
to the emission time of an optical phonon of about 150 fs, the carriers populate subsequent replicas at 
increasing times and, therefore, the efficiency of carrier-carrier scattering processes in broadening the 
peaks increases. However, the measured spectra exhibit a different behavior which is quantitatively 
reproduced by the SBE model: Already the unrelaxed peak exhibits a strong increase with increasing 
density and the density dependence of all replicas is approximately the same. The reason is that the 
broadening of the generation process as discussed |above obviously dominates over the broadening due 
to subsequent scattering processes. As discussed inL3, this is due to the fact that the generation process 
is strongly inffuenced by the decay of the interband polarization which, in turn, is due to electron plus 
hole scattering, while the broadening of the electron distribution during the relaxation process is due 
to electron scattering only. |^ 

The above theoretical and experimental analysis constitutes a clear demonstration of the importance 
of dephasing processes for the analysis of luminescence spectra. The density dependence of the spectra 
in the interesting transition region between carrier-phonon and carrier-carrier dominated dynamics 
can only be explained by a coherent modeling of the carrier generation including the dynamics of the 
interband polarization. The main mechanism determining the width of the peaks is the broadening 
of the generation rate. It can be shown that for band-to-acceptor spectea the broadening of the 
recombination process, neglected in the present model, plays a minor roleu. However, for band-to- 
band spectra also this phenomenon should be taken into account. For the case discussed above the 
background due to band-to-band luminescence has been found to be of negligible importance^. 



^ For the case of bulk GaAs considered here, due to the different electron and heavy-hole effective masses, the 
hole-hole scattering is about 5 times larger than the electron-electron one. 
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IV. BLOCH OSCILLATIONS AND WANNIER-STARK LOCALIZATION IN 

SUPERLATTICES 



Ever since the initial applications of quantum mechanics to the dynamics of electrons in solids, the 
analysis of Bloch electrons moving in a homogeneous-electric field has been of central importance. 
By employing semiclassical arguments, in 1928 BlochE3 demonstrated that, a wave packet given by 
a superposition of single-band states peaked about some quasimomentum, Kk, moves with a group 
velocity given by the gradient of the energy-band function with respect to the quasimomentum and 
that the rate of change of the quasimomentum is proportional to the applied field F. This is often 
referred to as the "acceleration theorem" : 

hk^eF . (56) 

Thus, in the absence of interband tunneling and scattering processes, the quasimomentum of a Bloch 
electron in a homogeneous and static electric field will be uniformly accelerated into the next Brillouin 
zone in a repeated-zone scheme (or equivalently undergoes an Umklapp process back in to the first 
zone). The corresponding motion of the Bloch electron through the periodic energy-band structure, 
shown in figure U, is called "Bloch oscillation"; It is characterized by an oscillation period tb — h/eFd, 
where d denotes the lattice periodicity in the field direction. 

There are two mechanisms impeding a fully periodic motion: interband tunneling and scattering pro- 
cesses. Interband tunneling is an intricate problem and still at the center of a continuing debate. Early 
calculations of the tunneling probability into other bapds in which the electric field is represented by a 
time-independent scalar potential were made by ZenerEJ using a Wentzel-Kramers-Brillouin generaliza- 
tion of Bkjch functions ,|-b.y HoustonCJ using accelerated Bloch states (Houston states), and subsequently 
by KaneEil and Argyrea23 who employed the crystal-momentum representation. Their calculations lead 
to the conclusion that the tunneling rate per Bloch period is much less than unity for electric fields up 
to lO^V/cm for typical band parameters corresponding to elemental or compound semiconductors. 

Despite the apparent agreement among these calculations, the validity of employing the crystal- 
momentum representation or Houston functions to describe electrons moving in a non-periodic (crystal 
plus external fielcU-potential has been disputed. The starting point of the controversy was the original 
paper by WannieiO. He pointed out that, due to the translational symmetry of the crystal potential, if 
(/)(r) is an eigenfunction of the scalar-potential Hamiltonian (corresponding to the perfect crystal plus 
the external field) with eigenvalue e, then any (jiir+nd) is also an eigenfunction with eigenvalue e + nAe, 
where Ae = eFd is the so-called Wannier-Stark splitting (d being the primitive lattice vector along the 
field direction). He concluded that the translational symmetry of the crystal gives rise to a discrete 
energy spectrum, the so-called Wannier-Stark ladder. The states corresponding to these equidistantly 
spaced levels are localized states, as schematically shown in figure || for the case of a semiconductor 
superlattice. |— ■ 

The existence of such energy quantization was disputed by Zakcj, who pointed out that for the 
case of an infinite crystal the scalar potential — F • r is not bounded, which implies a continuous energy 
spectrum. Thus, the main point of the controversy was related to the existence (or absence) of Wannier- 
Stark ladders. More pr(5cisely, the point was to decide if interband tunneling (neglected in the original 
calculation by WannieiEj) is so strong to destroy the Wannier-Stark energy quantization (and the 
corresponding Bloch oscillations) or not. 

It is only during the last decade that this controversy came to an end. From a theoretical point of view, 
most of the formal problems related to the non-periodic nature of the scalar potential (superimposed 
to the periodic, crystal potential) were finally removed by using a vector-potential representation of the 
applied fieldEJO. Within such vector-potential picture, upper boundaries for the interband tunneling 
probability have been established at a rigorous level, whieb-show that an electron may execute a number 
of Bloch oscillations before tunnelins-oiii. of the bandESEZi, in qualitatively good agreement with the 
earlier predictions of Zener and KancSo. 

The second mechanism impeding a fully periodic motion is scattering by phonons, impurities, etc. 
(see figure ^ . This results in lifetimes shorter than the Bloch period tb for all reasonable values of the 
electric field, so that Bloch oscillations should not be observable in conventional solids. 

In superlattices, however, the situation is much more favourable because of thejSinaller Bloch period 
Tb resulting from the small width of the mini-Brillouin zone in the field directiorO. 

Indeed, the existence of Wannier-Stark laddeifs as well as Bloch oscillations in superlattices has been 
confirmed by a number of recent experimental. The photoluminescence and photocurrent measure- 
ments of the biased GaAs/GaAlAs superlattices performed by Mendez and coworkersLj, together with 



the electroluminescence experiments by Voisin and coworkersE£il, provided the earliest evidence of 
field-induced Wannier-Stark ladders in superlattices. A few years later, Feldmann and coworkerst£il 
were able to measure Bloch oscillations in the rtiiiie domain through a four-wave-mixing experiment 
originally suggested by von Plessen and Thomasll23. A detailed analysis of the Bloch oscillations in the 
four- wav e , - m i ixin g signal (which reflects the interband dynamics) has been also performed by Leo and 



coworkergl23ll2: 

In addition to the above interband-polarization analysis, Bloch oscillations have been also detected 
by monitoring the intraband polarization which, in turn, is reflected by anisotropic changes in the 
refractive indexfJ. Measurements ^aa^d pn transmittive electrooptic sampling (TEOS) have been per- 



aasipd 

formed by Dekorsy and coworkergi2§E2^ 
direct detection of the TeraHertz (THz) 



Finally Bloch oscillations were recentljj- 
radiation in semiconductor superlatticesc 



asured through a 



A. Two equivalent pictures 



Let us now apply the theoretical approach presented in section || to the case of a semiconductor 
superlattice in the presence of an uniform (space-independent) electric field. The non-interacting car- 
riers within the superlattice crystal will then be described by the Hamiltonian H° in equation (|23|), 
where now the electrodynamic potentials A2 and (m the following simply denoted with A and 1^9) 
correspond to a homogeneo us el ectric field E2(r,t) — F(t). 

As pointed out in section II A , the natural quantum-mechanical representation is given by the eigen- 
states of this Hamiltonian: 



(-iW,-fA(r,t))^ 



-l-e(^(r,i) + y'(r) 



'/>n(r) = (r) • 



(57) 



However, due to the gauge freedom discussed in section [I A, there is an infinite number of possible 
combinations of A and Lp — and therefore of possible Hamiltonians — which describe the same homoge- 
neous electric field r(i). In particular, one can identify two independent choices: the vector-potential 
gauge 



A(r,i) = -c [ F{t')dt' , (p{r,t) 







and the scalar-potential gauge 



A(r,i)=0, (^(r,t) = -F(t)T 



(58) 



(59) 



As shown inE2a, the two independent choices correspond, respectively, to the well known Bloch- 
oscillation and Wannier-Stark pictures. They simply reflect two equivalent quantum- mechanical repre- 
sentations and, therefore, any physical phenomenon can be described in both pictures. 



, the eigenfunctipns (/)„ in { pT\ ) are the so- 
As discussed ir£2a, such time dependent 



More specifically, within the vector-potential pictJire 
called accelerated Bloch states (or Houston states)L3t 
representation constitutes the natural basis for the description of Bloch oscillations, i.e. it provides a 
rigorous quantum-mechanical derivation of the acceleration theorem (56), thus showing that this is not 
a simple semiclassical result. ^ Within such representation, Bloch oscillations are fully described by the 
diagonal terms of the intraband density matrix (|3|) (semiclassical distribution functions). Therefore, 
non-diagonal elements, describing phase-coherence between different Bloch states, do not contribute to 
the intraminiband dynamics. However, they are of crucial importance for the description of intermini- 
band dynamics, i.e. field-induced Zener tunneling, which in this Bloch-state representation originates 
from the time variation of our basis states (see equation (4^)). 

On the contrary, within tka scalar-potential picture (|5£ ) , the eigenfunctions (/>„ in (^7|) are the well- 
known Wannier-Stark statesEJ. Contrary to the previous Bloch picture, within such representation the 
intraminiband Bloch dynamics originates from a quantum interference between different Wannier-Stark 
states, thus involving non-diagonal elements of the intraband density matrix (pq). 



®The acceleration theorem (^6|) and the corresponding Bloch-oscillation dynamics are usually regarded as a 
semiclassical result compared to the Wannier-Stark picture. On the contrary, they correspond to two different 
fully quantum-mechanical pictures. 
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B. Some simulated experiments 



In this section, we will rjejdejSK-fecent simulated experiments of the ultrafast carrier dynaiaica in 



semiconductor super latticesli22 111. They are based on a generalized Monte Carlo solutio nEj LJ of 



the set of kinetic equations (so-called semiconductor Bloch equations) derived in section HE. In 



this case, the Bloch representation discussed in section [TV A| has been employed limiting the set of 
interband density-matrix elements in ( p9[ ) to the diagonal ones, i.e. i = j. In additiQn.. incoherent 
scattering processes have been treated within the usual Markov limit as discussed inoElllIj. Within 
such approximation schemCj-tiie explicit form of the kinetic equations (Q) coincides with the SBE ( p7| ) 
(obtained for the bulk case)t23, provided we replace the wavevector k with hi^, v being the superlattice 
miniband index. 

In the simulated experiments reviewed here the following superlattice model has been employed: The 
energy dispersion and the corresponding wavefunctions along the growth direction (fc||) are computed 
within the well known Kronig-Penney modelo, while for the in-plane direction (fcj^) an effective- mass 
model has been used. Starting from these three-dimensional wavefunctions 0^1/' various carrier- 
carrier as well as carrier-phonon matrix elements are numerically computed (see equations ( ^5| ) and 
(p7[)). They are, in general, functions of the various miniband indices and depend separately on fc|| and 
A;_L, thus fully reflecting the anisotropic nature of the superlattice structure. 

Only coupling to GaAs bulk phonons has been considered. This, of course, is a simplifying approxi- 
mation which neglects any superlattice effect on the phonon dispersion, such , as , confinement of optical 
modes in the wells and in the barriers, and the presence of interface modesEi£rEi3. However, while these 
modifications have important consequences for phonon spectroscopies (like Raman scattering) , they are 
far less decisive for transport phenomena. ^ 

We will start discussing the scattering-induced damping of Bloch oscillations. In particular, we-fp 
show that in the low-density limit this damping is mainly determined by optical-phonon scatterin,"'^*" 



while at high densities the main-mechanism responsible for the suppression of Bloch oscillations is found 



to be carrier-carrier scatteringtii . 



This Bloch-oscillation analysis in the ti me do main is also confirmed by its counterpart in the fre- 



quency domain. As pointed out in section IV A , the presence of Bloch oscillations, due to a negligible 
scattering dynamics, should correspond to Wannier-Stark energy quantization. This is confirmed by 
the simulatecUiBtical-absorption spectra, which clearly show the presence of the field-induced Wannier- 
Stark laddersEJ. 



1. Bloch-oscillation analysis 

All t|b£_simulated experiments presented in this section refer to the superlattice structure considered 
in Ref.Eill: 111 A GaAs wells and 17 A Alo.3Gao.7As barriers. Eor such a structure there has been 
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experimental evidence for a THz-emission from Bloch oscillationsli! 

In the first set of simulated experiments an initial distribution of photoexcited carriers (electron-hole 
pairs) is generated by a 100 fs Gaussian laser pulse in resonance with the first- miniband exciton {hujL ~ 
1540 meV). The strength of the applied electric field is assumed to be 4 kV/cm, which corresponds to 
a Bloch period tb = h/eFd of about 800 fs. 

In the low-density limit (corresponding to a weak laser excitation) , incoherent scattering processes do 
not alter the |Bl o ch rOscillation dynaiBics-Jhis is due to the following reasons: In agreement with recent 
experimentalt2ZHliZl and theoreticalll2a~Eiil investigations, for superlattices characterized by a miniband 
width smaller than the LO-phonon energy — as for the structure considered here — and for laser excita- 
tions close to the band gap, at low temperature carrier-phonon scattering is not permitted. Moreover, 
in this low-density regime carrier-carrier scattering plays no role: Due to the quasi-elastic nature of 
Coulomb collisions, in the low-density limit the majority of the scattering processes is characterized 
by a very small momentum transfer; As a consequence, the momentum relaxation along the growth 
direction is negligible. As a result, on this picosecond time-scale the carrier system exhibits a coherent 
Bloch-oscillation dynamics, i.e. a negligible scattering-induced dephasing. This can be clearly seen 



^Indeed, by now it is well knownElS that the total scattering rates are sufficiently well reproduced if the phonon 



spectrum is assumed to be bulk-like. 



from the time evolution of the carrier distribution as a function of k\\ (i.e. averaged over k±) shown 
in figure ^. During the laser photoexcitation (t = 0) the carriers are generated around k\\ = 0, where 
the transitions are close to resonance with the laser excitation. According to the acceleration theorem 
(^6|), the electrons are then shifted in k-space. When the carriers reach the border of the first Brillouin 
zone they are Bragg reflected. After about 800 fs, corresponding to the Bloch period tb, the carriers 
have completed one oscillation in k-space. As expected, the carriers execute Bloch oscillations without 
loosing the synchronism of their motion by scattering. This is again shown in figure ^, where we have 
plotted: (b) the mean kinetic energy, (c) the current, and (d) its time derivative which is proportional 
to the emitted far field, i.e. the THz-radiation. (It can be shown that, by neglecting Zener tunneling, 
the intraband polarization P'^/'' in equation (|5^) is proportional to the current.) All these three quan- 
tities exhibit oscillations characterized by the same Bloch period tb- Due to the finite width of the 
carrier distribution in k-space (see figure ^(a)), the amplitude of the oscillations of the kinetic energy is 
somewhat smaller than the miniband width. Since for this excitation condition the scattering-induced 
dephasing is negligible, the oscillations of the current are symmetric around zero, which implies that 
the time average of the current is equal to zero, i.e. no dissipation. 

As already pointed out, this ideal Bloch-oscillation regime is typical of a laser excitation close to gap 
in the low-density limit. Let us now discuss, still at low densities, the case of a laser photoexcitation 
high in the band. Figure ^a) shows the THzpsignal as obtained from a set of simulated experiments 
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corresponding to different laser excitationslliil. The different traces correspond to the emitted THz- 
signal for increasing excitation energies. We clearly notice the presence of Bloch oscillations in all 
cases. However, the oscillation amplitude and decay (effective damping) is excitation-dependent. 
For the case of a laser excitation resonant with the first-miniband cxciton considered above (see figure 
, we have a strong THz-signal. The amplitude of the signal decreases when the excitation energy is 
increased. Additionally, there are also some small changes in the phase of the oscillations, which are 
induced by the electron-LO phonon scattering. 

When the laser energy comes into resonance with the transitions between the second electron and 
hole minibands {hujL ~ 1625 meV), the amplitude of the THz-signal increases again. The correspond- 
ing THz-transients show an initial part, which is strongly damped and some oscillations for longer 
times that are much less damped. For a better understanding of these results, we show in figure 0(b) 
the individual THz-signals, originating from the two electron and two heavy-hole minibands for the 
excitation with Ticu — 1640 meV. The Bloch oscillations performed by the electrons within the jSetfpt^i 
miniband are strongly damped due to intra- and interminiband LO-phonon scattering proccsseai 



Since the width of this second miniband (45 meV) is somewhat larger than the LO-phonon energy, 
also intraminiband scattering is possible, whenever the electrons are accelerated into the high-energy 
region of the miniband. The THz-signal originating from electrons within the first miniband shows an 
oscillatory behavior with a small amplitude and a phase which is determined by the time the electrons 
need to relax down to the bottom of the band. 

At the same time, the holes in both minibands exhibit undamped Bloch oscillations, since the mini- 
bands are so close in energy that for these excitation conditions no LO-phonon emission can occur. 
The analysis shows that at early times the THz-signal is mainly determined by the electrons within 
the second miniband. At later times the observed signal is due to heavy holes and electrons within the 
first miniband. 

The above theoretical analysis closely resembles-experimental observations obtained for a superlattice 
structure very similar to the one modelled hereL£3. In these experiments, evidence for THz-emission 
from Bloch oscillations has been reported. For some excitation conditions these oscillations are associ- 
ated with resonant excitation of the second miniband. The general behavior of the magnitude of the 
signals, the oscillations and the damping are close to the results shown in figure |^. 

Finally, in order to study the density dependence of the Bloch-oscillation damping, let us go back to 
the case of laser excitations close to gap. Figure ||(a) shows the total (electrons plus holes) THz-radiation 
as a function of time for three different carrier densities. With increasing carrier density, carrier-carrier 
scattering becomes more and more important: Due to Coulomb screening, the momentum transfer in a 
carrier-carrier scattering increases (its typical value being comparable with the screening wavevector). 
This can be seen in figure ^(a), where for increasing carrier densities we realize an increasing damping 
of the THz-signal. However, also for the highest carrier density considered here we deal with a damping 
time of the order of 700 fs, which is much larger than the typical dephasing time, i.e. the decay time of 
the interband polarization, associated with carrier-carrier scattering. The dephasing time is typically 
investigated by means of .fnijr- wave-mixing (FWM) measurements and such multi-pulse experiments 
can be simulated as wellEilO. From a theoretical point of view, a qualitative estiiaiate of the dephasing 
time is given by the decay time of the "incoherently summed" polarization (ISP)e3. Figure ^(b) shows 
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such ISP as a function of time for the same three carrier densities of figure §(a). As expected, the decay 
times are always much smaller than the corresponding damping times of the THz-signfils (note the 



different time-scale in figures ^(a) and (b)). This difference, discussed in more detail inEJina, can be 
understood as follows: The fast decay times of figure ||(b) reflect the interband dephasing, i.e. the sum 
of the electron and hole scattering rates. In particular, for the Coulomb interaction this means the sum 
of electron-electron, electron- hole, and hole-hole scattering. As for the case of bulk GaAs discussed in 
section III, this last contribution is known to dominate and determines the dephasing time-scale. On 
the other hand, the total THz-radiation in figure ||(a) is the sum of the electron and hole contributions. 
However, due to the small value of the hole miniband width compared to the electron one, the electron 
contribution will dominate. This means that the THz damping in figure |^(a) mainly reflects the 
damping of the electron contribution. This decay, in turn, reflects the intraband dephasing of electrons 
which is due to electron-electron and electron-hole scattering only, i.e. no hole-hole contributions. 

From the above analysis we can conclude that the decay time of the THz-radiation due to carrier- 
carrier scattering differs considerably from the corresponding dephasing times obtained from a FWM 
experiment: The first one is a measurement of the intraband dephasing while the second one reflects 
the interband dephasing. |^ 



2. Optical-absorption analysis 



Let us now discuss the frequency-domain counterpart of the Bloch-oscillation picture considered so 
far. Similar to what happens in the time domain, for sufficiently high electric fields, i.e. when the Bloch 
period tb = h/eFd becomes smaller than the dephasing time, the optical spectra of the superlattice are 
expected to exhibit the frequency-domain coun terpart of the Bloch oscillations, i.e. the Wannier-Stark 



energy quantization discussed in section [IV A| . In the absence of Coulomb interaction, the Wannier- 
Stark ladder absorption increases as a function of the photon energy in a step-like fashion. These steps 
are equidistantly spaced. Their spacing, named Wannier-Stark splitting, is proportional to the applied 
electric field. 

The simulated linear-absorption spectra correspondingjto-|a superlattice structure with 95 A GaAs 
wells and 15 A Alo.3Gao.7As barriers are shown in figure ffij. As we can see, the Coulomb interaction 
gives rise to excitonic peaks in the absorption spectra and introduces couplings between these Wannier- 
Stark states. Such exciton peaksy, which are no longer equidistantly spaced, are often referred to as 



excitonic Wannier-Stark ladderslia of the superlattice. 



Since for the superlattice structure considered in this simulated experimentli£3iH3 the combined 
miniband width is larger than the typical two- and three-dimensional exciton binding energies, it is 
possible to investigate the quasi-three dimensional absorption behavior of the delocalized miniband 
states as well as localization effects induced by the electric field. 

For the free-field case, the electron and hole states are completely delocalized in our three-dimensional 
k-space. The perturbation induced by the application of a low field (here « 5 kV/cm), couples the 
states along the field direction and in the spectra the Franz-Keldysh effect, well known from bulk 
materialsEl, appears: one clearly notices oscillations which increase in amplitude with the field and shift 
with F'^/^ from the n = and n = 1 levels toward the center of the combined miniband. 

For increasing field the potential drop over the distance of a few quantum wells eventually exceeds the 
miniband width and the electronic states become more and more localized. Despite the field-induced 
energy difference neFd, the superlattice potential is equal for quantum wells separated by nd. Therefore, 
the spectra decouple into a series of peaks corresponding to the excitonic ground states of the individual 
electron- hole Wannier-Stark levels. Each Wannier-Stark transition contributes to the absorption with a 
pronounced 1 — s exciton peak, plus higher bound exciton and continuum states. The oscillator strength 
of a transition n is proportional to the overlap between electron and hole wavefunctions centered at 
quantum wells n' and n + n', respectively. The analysis shows that this oscillator strength is almost 
exclusively determined by the amplitude of the electron wavefunction in quantum well n' + n since 



®We stress that this difference between intraband and interband dephasing in superlattices is the same discussed 
in section III for the case of bulk semiconductors, where the broadening of the photoexcited carrier distribution 



is mainly determined by the decay of the interband polarization (interband dephasing) while the subsequent 
energy broadening of the electron distribution is due to electron scattering only (intraband dephasing) (see 
figure 1). 



on 



for fields in the Wannier-Stark regime the hole wavefunctions are almost completely localized over one 
quantmn well due to their high effective mass (see figure ||). Thus, the oscillator strengths of transitions 
to higher |n| become smaller with increasing \n\ and field. 

At high fields (here >« 8 kV/cm) the separation between the peaks is almost equal to neFd. For 
example, the peak of the n = transition which is shifted by the Wannier-Stark exciton binding energy 
with respect to the center of the combined miniband, demonstrates that the increasing localization also 
increases the exciton binding energy. This increased excitonic binding reflects the gradual transition 



from a three- to a two-dimensional behavior, discussed in detail inc 

For intermediate fields there is an interplay between the Wannier-Stark and the Franz-Keldysh effect. 
Coming from high fields, first the Wannier-Stark peaks are modulated by the Franz-Keldysh oscillations. 
However, as soon as the separation eFd between neighboring peaks becomes smaller than their spectral 
widths, the peaks can no longer be resolved individually so that only the Franz-Keldysh structure 
remains. 



V. COULOMB-CORRELATION EFFECTS IN SEMICONDUCTOR QUANTUM WIRES 



The importance of Coulomb-correlation effects in opticqiL-^pectra of semiconductors and their de- 
More recently, increasing interest has 



pendence on dimensionality has now been long-cecognizedti^ . 

been devoted to one-dimensional (ID) systemsE£3, prompted by promising advances in quantum-wire 
fabrication and application, e.g. quantum-wire lasers. The main goal is to achieve structures with 
improved optical efficiency as compared to two-dimensional (2D) and three-dimensional (3D) ones. A 
common argument in favour of this effort is based on the well known van Hove divergence in the ID 
joint density-of-states (DOS), which is expected to give rise to very sharp peaks in the optical spectra 
of ID structures. Such prediction is however based on free-carrier properties of ideal ID systems and 

it ignores Coulomb-correlation effects. .j—. 

Early theoretical investigations by Ogawa and coworkersE£llE£3, based on a ID single-subband model 
of the wires, showed that the inverse-square-root singularity in the ID DOS at the band edge is smoothed 
when excitonic effects are taken into account. 



In this section, we will review recent theoretical resultsEH5ll23 based on a full 3D description of 
realistic quantum- wire structures made available by the present tecknology, such as structures obtained 
by epitaxial growth on non-planar substrates (V-shaped wires )EH3 or by cleaved-edge quantum well 
overgrowth (T-shaped wires)tH3. On the one hand, this approach allows an accurate determination of 
ID-exciton binding energies; on the other hand, it allows us to investigate also the non linear (gain) 
regime. In both cases, one finds a strong suppressiati, pl-the ID band-edge singularity, in agreement 
with previous results based , on ,, sin iiplified ID modeli 
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This theoretical approachllH3llH3 is based on the general kinetic theory presented in section |l[ which 



allows a full 3D description of Coulomb interaction within a multiminiband scheme. In particular, we 
focus on the quasi-equilibrium regime where the solution of the coupled kinetic equations (Q) simply 
reduces to the solution of the interband-polarization equatipii. This is performed by direct numerical 
evaluation of the polarization eigenvalues and eigenvectorsEHJ, which fully determine the absorption 
spectrum as well as the exciton wavefunctions. The main ingredients entering this calculation are 
the single-particle energies and-wavefunctions, obtained numerically for the 2D confinement potential 



deduced, e.g., from TEM as intHS 



The above theoretical scheme has been applied to realistic V- and T-sha]:fed| wire structures. In 

whose cross-section 



particular, here we show results for the GaAs/AlGaAs V-wire structure of Ref.P^ 
is shown in figure UU. 



A. Linear response: excitonic-absorption regime 



Let us start considering the optical response of the system in the low-density limit. In figure |Tl| we 
show the linear-absorption spectra obtained when taking into account the lowest wire transition only. 
Results of our Coulomb-correlated (CC) approach are compared to those of the free-carrier (FC) model . 
As we can see, electron-hole correlation introduces two important effects: First, the excitonic peak arises 
below the onset of jtlie continuum, with a binding energy of about 12 meV, in excellent agreement with 
recent experimentgi23. Second, one finds a strong-si 
analysis of the physical origin of such suppression 



ession of the ID DOS singularity. A detailed 
EHJ has shown that the quantity which is mainly 
modified by CC is the oscillator strength (OS). In flgure 12^(a) the ratio between the CC and FC OS is 
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plotted as a function of excess energy (solid line). This ratio is always less than one and goes to zero 
at the band edge, and reflects a sort of hole in the electron-hole correlation function g{z) as shown in 
figure ^(b). Such vanishing behaviour is found to dominate the ID DOS singularity and, as a result, 
the absorption spectrum exhibits a regular behaviour at the band edge (solid line in figure pi] ) . 

The above analysis shows that also for realistic quantum-wire structures electron-hole correlation 
leads to a strong suppression of the ID band-edge singularity in the linear-absorption spectrum, contrary 
to the 2D and 3D cases. 



B. Non-linear response: gain regime 



Most of the potential quastum-wire applications, i.e. ID lasers and modulators, operate in strongly 



non-linear-response regimealHj. In such conditions, the above linear-response analysis has to be gen- 
eralized taking into account additional factors as: (i) screening effects, (ii) band renormalization, (iii) 
space-phase filling. 

Figure reports quantitative results for non-linear absorption spectra of realistic V-shaped wire 
structures at different carrier densities at room temperature. As a reference we also show the results 
obtained by including the lowest subband only [figure [T^(a)]. In the low-density limit (case A: n — 10^ 
cm~^) we clearly recognize the exciton peak. With increasing carrier density, the strength of the 
excitonic absorption decreases due to phase-space filling and screening of the attractive electron-hole 
interaction, and moreover the band renormalization leads to a red-shift of the continuum. At a density 
of 4 * 10^ cm~^ (case D) the spectrum already exhibits a negative region corresponding to stimulated 
emission, i.e. gain regime. As desired, the well pronounced gain spectrum extends over a limited energy 
region (smaller than the thermal energy); However, its shape differs considerably from the ideal FC one 
(curve marked with diamonds in the same figure). In particular, the band-edge singularity in the ideal 
FC gain spectrum is clearly smeared out by electron-hole correlation. The overall effect is a broader 
and less pronounced gain region. 

Finally, figure ^(b) shows the non-linear spectra corresponding to the realistic case of a 12-subband 
V-shaped wire. In comparison with the single-subband case [figure |l^(a)] , the multisubband nature 
is found to play an important role in modifying the typical shape of the gain spectra, which for both 
CC and FC models result to extend over a range much larger than that of the single-subband case 
for the present wire geometry. In addition, the Coulomb-induced suppression of the single-subband 
singularities, here also due to intersubband-coupling effects, tends to reduce the residual structures 
in the gain profile. Therefore, even in the ideal case of a quantum wire with negligible disorder and 
scattering-induced broadening, our analysis indicates that, for the typical structure considered, the 
shape of the absorption spectra over the whole density range differs considerably from the sharp FC 
spectrum of figure |l^. 

This tells us that, in order to obtain sharp gain profiles, one of the basic steps in quantum- wire tech- 
nology is to produce structures with increased subband splitting. However, the disorder-induced inho- 
mogeneous broadening, not considered here, is known to increase significantly the spectral broadening 
and this effect is expected to increase with increasing subband splitting. Therefore, extremely-high- 
quality structures (e.g. single-monolayer control) seem to be the only possible candidates for successful 
quantum-wire applications. 



VI. SUMMARY AND CONCLUSIONS 



A review of coherent phenomena in optically-excited semiconductors has been presented. Our analysis 
has allowed the identification of two classes of phenomena: optically-induced and Coulomb-induced 
coherent phenomena. We have shown that both classes can be described in terms of a unique theoretical 
framework based on the density-matrix formalism. Due to its generality, such quantum-kinetic approach 
allows a realistic description of coherent as well as incoherent, i.e. phase-breaking, processes, thus 
providing quantitative information on the coupled — coherent vs. incoherent — carrier dynamics in 
photoexcited semiconductors. 

The primary goal of the paper was to discuss the concept of quantum-mechanical phase coherence 
as well as its relevance and implications on semiconductor physics and technology. In particular, 
we have analyzed the dominant role played by optically induced phase coherence on the process of 
ultrafast carrier photogeneration. We have then discussed typical field-induced coherent phenomena in 
semiconductor superlattices, e.g. Bloch oscillations and the corresponding THz radiation, as well as 
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their dephasing dynamics. Finally, we have analyzed the dominant role played by Coulomb correlation 
on the linear and non-linear optical spectra of realistic quantum- wire structures, namely the Coulomb- 
induced suppression of ideal ID-band-cdgc singularities. 

Our analysis shows that the conventional separation between coherent and incoherent regimes is 
no longer valid for most of the recent ultrafast optical experiments in semiconductors. The reason 
is that on such extremely short timc-scalcs the carrier dynamics is the result of a strong interplay 
between coherence and relaxation, thus preventing a time-scale separation between photogeneration 
and relaxation dynamics. Similar considerations apply to the various Coulomb-induced phenomena: 
The strong modifications induced by Coulomb correlation on the optical response of low- dimensional 
semiconductors play a significant role on the phase-breaking relaxation dynamics as well. It is thus 
evident that any theoretical analysis of the optical response of semiconductor heterostructures must 
provide a proper description of both coherent and incoherent phenomena as well as of their mutual 
coupling on the same kinetic level. 
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FIG. 1. Generation rate of electrons as a function of k at different times and densities. After Ref0. 

FIG. 2. Galculated and measured band-to-acceptor spectra for different carrier densities. After RefEl. 

FIG. 3. Full width at half maximum of the unrelaxed peak and the first and second phonon replica as 
a function of carrier density obtained from the calculated iMnd measured spectra. The lines connecting the 
theoretical values are meant as guides to the eye. After RefEj. 

FIG. 4. Schematic illustration of the field-induced coherent motion of an electronic wavepacket initially 
created at the bottom of a miniband. Here, the wjidljh of the miniband exceeds the LO-phonon energy E^o, so 
that LO-phonon scattering is possible. After Ref. 



FIG. 5. Schematic representation of the transitioipifrom the valence to the conduction band of a superlattice 
in the Wannier-Stark localization regime. After Ref.tj. 

FIG. 6. Full Bloch-oscillation dynamics corresponding to a laser photoexcitation resonant with the 
first-miniband exciton. (a) Time evolution of the electron distribution as a function of fey . (b) Average ki- 
netic energy, (c) current, and (d) THz-signal corresponding to the Bloch oscillations in (a). 



FIG. 7. (a) Total THz-signals for eight different spectral positions of the exciting laser pulse: 
1540, 1560, . . . , 1680 meV (from bottom to top), (b) Individual THz-signal of the eleqteens and holes in the 
different bands for a central spectral position of the laser pulse of 1640 meV. After Ref. 
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FIG. 8. (a) To1jaJ-|THz-radiation as a function of time; (b) Incoherently-summed polarization as a function 
of time. After ReLE^. 

FIG. 9. Absorption spectra for various static applied electric fields for a GaAs/Alo.sGao.rAs superlattice 
(well (barrier) width 95 (15) A). The vertical displacements between any two spectra is proportional to the 
difference of the corresponding fields. The Wannier-Stapietransitions are labeled by numbers, the lower (higher) 
edge of the combined miniband by Eq {E\). After Ref. 



FIG. 10. Typical cross-section ol V-grooved wires derived Irom TEM micrographs and used to define the 
confinement potential entering the 2D single-particle Schrodinger equation. The wire width along y at the apex 
of the V is about 10 nm. The dashed lines identifiy the region (about 70 x 25 nm), where single-particle and 
excitonic wavefunctions are localized. After RefJlHj. 

FIG. 11. Linear-absorption spectra of realistic V-shaped quantum-wire structures obtained by including 
the first electron and h|ei« subbands only. Solid line: Coulomb correlated (CC) result; Dashed line: free-carrier 
(PC) result. After Ref.H. 



FIG. 12. (a) OS ratio and DOS vs. excess energy; (b) electron-hole correlation function. After Ref. 



FIGi 13. Non-linear absorption spectra: (a) single-subband case; (b) realistic 12-subband case. After 
Ref.E!l. 
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